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Abstract 

A dynamic model for a random network evolving in continuous time is 
defined where new vertices are born and existing vertices may die. The 
fitness of a vertex is defined as the accumulated in-degree of the ver¬ 
tex and a new vertex is connected to an existing vertex with probability 
proportional to a function b of the fitness of the existing vertex. Further¬ 
more, a vertex dies at a rate given by a function d of its fitness. Using 
results from the theory of general branching processes, an expression for 
the asymptotic empirical fitness distribution {pk} is derived and analyzed 
for a number of specific choices of b and d. When b{i) = i+a and d{i) = fi 
- that is, linear preferential attachment for the newborn and random 
deaths - then pk ^ When b{i) = i + \ and d(i) = j3{i + 1), with 

/3 < 1, then pk ^ {1+ /3)“^, that is, if also the death rate is proportional 
to the fitness, then the power law distribution is lost. Furthermore, when 
b{i) = i 1 and d{i) = P{i + !)'>', with /3,7 < 1, then logpfe - 

a stretched exponential distribution. The momentaneous in-degrees are 
also studied and simulations suggest that their behaviour is qualitatively 
similar to that of the fitnesses. 

Keywords: Branching process, random network, preferential attachment, 
degree distribution, power law distribution. 
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1 Introduction 

Empirical studies on real networks have revealed that many of them exhibit 
features that are not captured by the classical Erdos-Renyi graph. In par¬ 
ticular, many networks tend to have a quite heavy tailed degree distribution, 
often described by a power law, that is, the fraction of vertices with degree k 
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decays as for some exponent r. To capture this, a number of new graph 
models have been developed that allow for this type of degree distribution; see 
for instance the seminal paper by Bollobas et al. [g (the setup described there 
includes many previously studied models as special cases) and the references 
therein. 

A model type that has received a lot of attention is based on preferen¬ 
tial attachment: new vertices are successively added to the network and are 
connected to existing vertices proportionally to a function f{k) = of the 
degree. The mechanism was introduced in the context of network modeling 
by Barabasi and Albert [3], who observed that, for 7 = 1 it seemed to lead to 
power law distributions with exponent r = 3 for the degrees. This was proved 
rigorously in [5]. Krapivsky and Redner m later derived heuristic results 
indicating that the regimes 7 < 1 and 7 > 1 are qualitatively different: for 
7 > 1 the degree distribution is degenerate in the sense that there is a single 
dominant vertex that is connected to almost every other vertex. This was 
confirmed in m- Recently Dereich and Morters [in] has studied the degree 
evolution in the sublinear regime 7 < 1 . 

Most preferential attachment models in the literature are defined in dis¬ 
crete time. Rudas et al. [20] however studies the model from a new perspec¬ 
tive, embedding it in continuous time and interpreting it as a continuous time 
branching process. This opens up for the use of well established results from 
the theory of branching processes. Compared to the traditional techniques 
for analyzing models with linear preferential attachment, the approach has 
the advantage that it applies for any 7 < 1. We mention also the paper by 
Athreya et al. | 2 ], where the linear case 7 = 1 is studied by aid of a different 
type of branching process embedding. 

The purpose of the current paper is to define and analyze a preferential 
attachment model embedded in continuous time where vertices are not just 
added but may also die, and where the death rate of a vertex is taken to be a 
function of its previous success, quantified by its accumulated degree. Using 
results from the theory of general branching processes, a characterization of the 
asymptotic distribution of the accumulated degrees is derived and analyzed 
in various special cases. The momentaneous degrees are also discussed and 
analyzed by means of simulations. This generalizes results from [20] . More 
specifically, the generalization compared to [ 20 | lies in that we work with a 
birth and death process instead of a pure birth process. We mention also the 
paper [1] by Britton and Lindholm, where a model is treated that is similar 
in spirit but not based on preferential attachment. 

Preferential attachment models incorporating deaths of vertices are much 
less studied than pure growth models, while real networks may well be affected 
by vertex emigration. Existing rigorous results are restricted to the case when 
vertices die randomly, that is, the death mechanism does not depend on fea¬ 
tures of the vertices (age, degree etc); see [7] and | 8 |. See also [16] and [9] for 
heuristic results, the latter concerning a situation where vertices die with a 
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probability inversely proportional to their degree. Concentration results are 
so far comparably weak, see e.g. O Lemma 4.1], and are lacking for models 
with more complicated death mechanisms. 

1.1 Definition of the model 

We now proceed to define the model. As already described, a vertex population 
evolving in continuous time is considered where existing vertices give birth to 
new vertices but may also die. More specifically, a vertex that is alive at time 
t and that has given birth to i new vertices during its life time, gives birth 
to an f + list vertex at rate b{i), that is, the probability that the vertex gives 
birth in the time interval {t,t + e) is 6(f)e + o(e). Furthermore, the vertex dies 
at rate d{i), that is, its probability of dying during (f, f + e) is d{i)e + o{e). Let 
^ be a random point process on [0, oo) with points (representing the births) 
generated according to the above rules. Starting from a single root vertex at 
time 0, each vertex x reproduces according to an independent copy of 
This means that the vertex population evolves as a birth and death process 
in continuous time and hence fits into the framework of general branching 
processes; see e.g. Jagers m Chapter 6]. Below we shall formulate conditions 
that guarantee that the process is supercritical, so that the vertex population 
indeed grows to infinity with positive probability. 

Write Zt for the set of vertices that are alive at time t, and, for x € Zt, let 
Af denote the total number of children that vertex x has produced up to time 
t (including children that are no longer alive), also referred to as the fitness 
of the vertex. The time until the next birth in the population is exponentially 
distributed with parameter Ylx^Zt ^ the time until the next death is 

exponentially distributed with parameter Ylx&Zt ^ (^f )• Hence the probability 
that the next event that occurs in the population is a birth is given by 

and the probability that the birth occurs at vertex x is proportional to b{Af). 
We shall see later that, under certain assumption on b and d, the ratio iJTj) 
converges almost surely to a number in (0,1). 

The vertex population described above naturally gives rise to a network 
structure: directed edges are created between a vertex and its children, point¬ 
ing from the child to the mother. When a vertex dies it turns into a ghost and 
all edges associated with it are turned into ghost links. According to the dy¬ 
namics, the rate at which a vertex produces/attracts new vertices depends on 
its previous success, measured by its in-degree, including ghost links. Weather 
or not it is realistic to include the ghost links depends on the setting. In many 
situations it is indeed natural to assume that the ability of a vertex to ac¬ 
quire new edges is affected by the total historical success of the vertex rather 
than just the contributions generated by vertices that are still present in the 
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network (think e.g. of coauthor-ship networks, sexual networks, recruitment 
networks etc). A model where ghost links are excluded in the definition of the 
fitness can be analyzed by the same means as the present model. The calcu¬ 
lations however become cumbersome and do not lead to explicit expressions. 
We elaborate on this and present simulation results in Section 5. 

1.2 Results 

Let p{t) be the density of that is, 

p{t) = lim e“^P((t, t + s) contains a point from 
£^0 

and dehne the Laplace transform 

poo 

p(X) = / e~^^p{t)dt. (2) 

Jo 

Clearly p{X) increases as A —0, and we write 

A = inf{A > 0 : p{\) < oo}. (3) 

Our first result is an expression for /5(A). 

Proposition 1.1 We have that 

OO k u/ ’\ 

Throughout, we shall assume that the functions b and d are such that 

p{\*) = 1 has a solution A* € (0,oo). (C.l) 

In fact, we shall assume that 

A* > A. (C.2) 

It is well-known from the theory of branching processes that, when A < oo, the 

condition (C.l) is equivalent to the process being supercritical. To be more 
precise, write ^{t) for the ^-measure of [0,t] and let Zt = \Zt\ be the number 
of vertices alive at time t. Then ¥{Zt —>■ oo) > 0 if and only if E[^(oo)] > 1 if 
and only if A* > 0; see m Chapter 6]. Writing 

p{X) = i°g(i-A+6w+d(i)) (4) 

k=0 

and using the Taylor expansion log(l—x) = —x+0(x^) reveals that a necessary 
condition for p(X) < oo is that 

A -|- d(i) _ 

^ A + 6(i) + d(i) ■ 
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and hence a necessary condition for A < oo (which is indeed necessary for 
(C.l)) is that 

1 + dii) 

^ bH) + dii) - 


Furthermore, a sufficient condition for having /5(A) < oo for all A > 0 (so that 
A = 0) is that 


1 + d{i) _^ 

— -—T > Cl ' tor some 7 < 1 . 

d{i) + b{i) - 


(5) 


Indeed, a bound of the same order then applies for and hence, by 

Taylor expansion and integral approximation, the exponent in (|3]) decays at 
least like making the sum convergent. Also note that (C.l) - or, more 

specihcally, the fact that A < 00 - guarantees that the process does not explode 
in finite time: a sufficient condition for non-explosion is that E[^(t)] < 00 for 
some f > 0 , see m Theorem 6.2.2], and, since p{t)dt = E[^((it)], it is clear 
that A cannot be finite if this does not hold. 

We are now ready to formulate the main theorem, which is a characteri¬ 
zation of the asymptotic htness distribution in the vertex population. To this 
end, recall that Zt is the number of alive vertices at time t and write for 
the number of alive vertices at time t with fitness k. Furthermore, let L be a 
random variable distributed as the life time of a vertex. 


Theorem 1.1 Assume thatb and d satisfy conditions (C.l) and (C.2). Then, 
for all k, on the event {Zt —>■ 00 }, almost surely 

1™ A=_^_pf_ m _ 

t^oo Zt A* -k b{k) + d{k) A* -k h{i) -k d{i) 

where C~^ = e“'’^**P(L > t)dt and the empty product arising for k = 0 is 

defined as 1. 

The rest of the paper is organized as follows. In Section 2 we state the 
results we shall need from the theory of branching process and, in Section 3, 
Proposition 11.11 and Theorem 11.11 are proved. Then, in Section 4, the limiting 
distribution in Theorem 1 1.1 1 is analyzed in more detail for a number of specific 
choices of h{i) and d{i). In Section 5, a formula is given for the the asymp¬ 
totic degree distribution in a network where the birth and death rates are 
based on the momentaneous in-degrees rather than the fitnesses (that is, the 
accumlated in-degrees). Unfortunately this expression is hard to analyze, but 
simulation results are presented that indicate that the qualitative behaviour 
of the distribution is the same as for the fitnesses. Finally, some possibilities 
for further work are mentioned in Section 6 . 
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2 Preliminaries on branching processes 


The corpus of literature on branching processes is vast and we do not intend to 
give a survey; for this, see Jagers m or the more expository book by Haccou 
et al. m- We shall however briefly explain the concepts that are needed 
to establish Theorem 0 These revolve around general branching processes 
counted with random characteristics. In a general branching process, the 
individuals reproduce independently of each other according to a random point 
process ^ on [0, 00 ) during a life-time L, where there are no restrictions on the 
dependence structure between ^ and L. An individual x can be represented 
as X = (ji,... ,jn) indicating that x is the jn'-th. child of the j„_i:th child ... 
of the ji:th child of the ancestor. Hence the space of possible individuals is 

00 

II = IJn”. 

i=0 

Write ax for the time when the individual x is born, that is, (Tq = 0 and, if x 
is the z:th child of the individual x', then ax = ax' + inf{t : > i}- 

Roughly, when counting a branching processes with a random characteris¬ 
tic, a random score is assigned to each individual and the quantity of interest 
is then the total score in the population. More specifically, let (f){t) be a 
real-valued stochastic process with = 0 for t < 0 giving the score of an 
individual of age t. The score of the individual x at time t is then given by 
4>x(t) = (j){t — ax) and the total score in the population at time t becomes 


Examples of common characteristics are (/)(t) = l{t < 0}, in which case zf 
is the number of individuals born up to time t, and 4>{t) = 1{0 < t < L}, 
in which case zf is the number of individuals that are alive at time t - this 
quantity we shall use the notation Zf for. These characteristics are both 
individual, meaning that ^x depends only on the process of descendants of 
X. Characteristics that involve information also about preceding generations 
are called relational. An example of such a characteristic that we shall be 
interested in below is when an individual is assigned the value 1 if its mother 
is still alive and 0 otherwise. 

The Laplace transform p(A) of the density p{t) of ^ and the related quantity 
A are defined as in ([2|) and ([3]). The solution A* to the equation /5(A) = 1 is 
called the Malthusian parameter and the process is subcritical, critical or 
supercritical depending on weather A* < 0, = 0 or > 0. A more frequent 
formulation of the definition of A* in the literature is as the solution to the 
equation 

poo 

/ e-^%adt)] = 1 , 

Jo 
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where, as in Section 1.1, ^{t) = C([0!^])- Since E[^(dt)] = p{t)dt, the formula¬ 
tions are clearly equivalent. 

Now, Theorem 11.11 will turn out to be a direct consequence of the following 
result concerning almost sure convergence of ratios of random characteristics. 

Theorem 2.1 (Nerman (1981)) Consider a supercritical branching process 
with Malthusian parameter X* > X and let cf) and ip be two individual charac¬ 
teristics satisfying 

E[sup{e“^*(/)(t)}] < oo for some X < X* (6) 

t 

and likewise for f:. Then, on the event {Zt —^ oo}, almost surely 

Zf . 0(A*) 

Zf ^(A*) 

where ^{X) = e“'’'*E[(/)(t)]dt. 

Remark 1. The characteristics we shall work with will in general be bounded 
and, for such characteristics, (l6|) is trivially satisfied. 

Remark 2. As pointed out in m pp 163], the conclusion of the theorem 
holds true also for relational characteristics provided that these depend only 
on a finite number of preceding generations; see [H] for details. 


3 Proofs 


We now return to the specific setting defined in Section 1.1, where ^ is the 
point process generated by the described birth and death rules. 

Proof of Proposition II.IL The expression for /5(A) follows from straightfor¬ 
ward calculations. First note that 

OO 

p{X) = V b{k) / e-^*P(^(t) = knL> t)dt, 

k=o Jo 

where trivially 


P(<f(t) = knL>t)= P(C(t) = k)- P(<f(t) = knL<t). 

Consider a generic vertex (the root vertex for instance) and write Tj for the 
time, counting from the birth of the vertex, when its z:th child is born, with 
Ti := oo if ^(oo) < i. Also, when Ti < oo, define = Tj+i — Ti and let 
be the time, counting from Tj, until the vertex dies. For A: = 0, we obtain 

roo poo 

/ e"^*P(^(t) = 0 n L > t)dt = / 

Jo Jo 

1 

A -|- 5(0) -|- d{0) 


7 






where e (^(o)+'^(o))* ig the probability that nothing happens in the process up 
to time t. To deal with k > 1, note that 

F{C{t) = k)= P(C(t) = k\Tk < oo)P(Tfc < oo), 

with 


= k\Tk < oo) = P(C(t) > k\Tk < oo) 

-F{^{t) > k + l\Tk+i < oo)P(Tfc+i < oolTfc < oo). 

Here 

F{m > k\Tk < CX)) = P < t\Tk < oo^ , 

where, conditionally on Tk < oo, the variables Ag,..., A^_^ are independent 
and A^ is exponentially distributed with parameter b{i) + d{i). Similarly, 

P(.^(t) = kn L < t) = P(C(i) = k n L < t\Tk < oo n Tk+i = oo) • 

F{Tk < oo n Tfc+i = oo) 


with 


^k-l 


P(.^(t) = knL < t\Tk < oonTfc_|_i = oo) = P E A^i+Af<t\Tk < cxnTfc+i 


,i=0 


where, conditionally on < oo fl Tk+i = oo, the variables Aj,..., A^_^, A 
are independent, A^ exponentially distributed with parameter b{i) + d{i) and 
A^ exponentially distributed with parameter b{k) + d{k). Also note that 


fc-i 


F{Tk < oo) = 


b{i) 




=: Qk 


F{Tk+i < oo\Tk < oo) = 


bik) 


P(Tfc < oo n Tk+i = oo) = 


d{k) 


b{k) + d{k) 

k-l 

n 


=: Qk 
b{i) 


b{k) + d{k) b{i) + d{i) 


w =: Qk- 


For independent exponential variables {Xi}, where Aj has parameter b{i) + 
d{i), write 


Ik = 


roo 

/ e-^¥(Ao + ... + Xk-i < t)dt. 

Jo 


Combining all of the above, we get that 

6 ( 0 ) 


/5(A) = 


A + 6(0) + d{0) 


+ 


'^b{k)(^qklk-{qkq'k + qk)Ik+iy (7) 


k=l 


?r Si- 



and all that remains is to note that 


1 r°° 

h = ^ fxo+...+Xk_i{t)dt 

_ 1 TT + d{i) 

A A + b{i) + d{i) ’ 

where the first equality follows from partial integration and the second from 
independence. Substituting this, and the expressions for and q'l, into 

(I?]) gives the desired formula for p(A). □ 

We proceed to prove the main result, Theorem 11.11 which, as already 
mentioned, follows from Theorem 12.11 

Proof of Theorem II.IL Consider the characteristic 

= /c n 0 < t < L}, 

that is, zf = is the number of alive vertices with k children at time t, 
and let V'(t) = 1{0 < t < L} so that zf = Zt is the number of alive vertices 
at time t. These characteristics are both bounded, and hence, since b and d 
are assumed to be such that (C.l) and {C.2) are satisfied, all assumptions of 
Theorem O are fulfilled. It follows that, on the event Zt —>■ oo, the ratio 
Z^jZt converges almost surely to the limit identified in Theorem 12.11 The 
expression in the nominator becomes 

l-QO 

0 (A*) = / = kf^L> t)dt, 

Jo 

and identical calculations as in the proof of Proposition II.II yields the formula 
in Theorem ll.il □ 


4 Examples 


We now investigate the behavior of the limiting fitness distribution in Theorem 
11.11 for a few specific choices of birth and death rates b and d. The limiting 
distribution is henceforth denoted by {pk}- Furthermore, for two sequences 
{ofc} and {bk} of real numbers, we write ~ bk to denote that ak/bk converges 
to a strictly positive constant as A: —>■ oo. 

Example 1. First consider the case when b{i) = i -\- a and d{i) = j3 with 
a > 0, that is, the birth rate is a linear function of the fitness and the vertices 
die at a constant rate. This means that a new vertex is connected to an 
existing vertex with a probability that depends linearly on the fitness of the 
existing vertex, and the vertices die at random. Using the relation 


k 


n(i+ c) 

i=0 


F(A: + l + c) 
F(c) 
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we get that 


„, . -p-r i -\- a r(A + q; + /3) r(A; + 1 + a) 

A + i + q: +/3 ^(q:) j-~^ V{k-\-\-\-ot-\-j3-\-X) 

k=0 *=u k=0 

Van der Hofstad et al. m Lemma 5.4] derives the formula 

^^r(A: + a) 1 /r(n + l + a) r(a) A 
^r(A: + c) 1 + a —cly r(n + c) r(c—1)/ 

and, by Stirlings formula, r(A; + c)/T{k) ~ as /c —oo. It follows that 




oo if A < 1 — /3; 
otherwise. 


Hence A* = 1 + a — /3, implying that the process is supercritical if /3 < a + 1. 
For the fitness distribution we get that 


Pk 


C 

k -\- X -\- 2(u 


n 


i-\- a 

i “h 1 “h 2 qi 


c 


r(A: + a) 

T{k + 2 + 2a) 


;,-( 2 +a)_ 


The limiting behavior of the fitness distribution is hence unaffected by the 
vertex death. This is explained by the fact that the vertices die randomly, 
that is, given that the next thing to occur is a vertex death, the vertex to 
die is chosen uniformly at random, and such a mechanism does not affect the 
tail behavior of the power law distribution generated by the birth mechanism. 
With a = 1 the expression for pk is in agreement with previous results for 
pure linear preferential attachment. 

Recall that, given the development of the process up to time t, the proba¬ 
bility that the next thing occurring is a vertex birth is given by m, where Af is 
the fitness of vertex x at time t. If the characteristics given by 4>x{t) = ) 

and ipxit) = d{Af_^^) satisfy the condition Q, then it follows from Theorem 
12.11 that the probability converges almost surely to a number in (0,1). This 
number is in general difficult to determine, but in the current case, denoting 
the total fitness Al°^ = YlxeZt we get 

Exez^HAf) _ Aj-^ + aZt 
Exez, ib (Af) + d (Af)) Af^ + (a + • 

The total fitness in the alive population can be obtained by counting the pop¬ 
ulation by the characteristic 'iA(t) = l{t > 0 and the mother of the individual 
is alive}. Hence, by Theorem 12.11 and the second remark following it, the ratio 
A\°^ jZt converges to 1 almost surely and it follows that 

1 -|- cr 

lim Pfthe next event is a vertex birth) =--. 

t^oo ^ 1 + a -b /3 
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This probability is strictly larger than 1/2 precisely when (5 < a + 1, that is, 
when the process is supercritical. 

Example 2. Next consider the case when b{i) = i + 1 and d{i) = /3(z + 1), 
that is, both the birth and the death rate are linear functions of the htness. 
In this case, by ([S|), we have p{X) < oo for all A > 0, and it is given by 


/5(A) 


i + 1 


OO k 

SS^(l + «i + (l+/9) 


E 


r(A: + 2 ) 


\ + l + ^ 

1 + ^ ^ r (A: + 1 + (1 +/?)' 

A + 1 + /? 


F 1 1 

’ ’ 1 + /3 ’1 + /3 


- 1 , 


where F denotes Gauss hypergeometric function; see Abramowitz and Stegun 
[U Chapter 15]. The equation p{\*) = 1 cannot be solved explicitly but a 
numerical solution is easily obtained for a given value of f5. The process is 
supercritical precisely when /3 < 1 , since the expected number of children 
produced by a vertex is given by 


OO k—1 


IE[C(oo)] = ^P(C(oo) > A:) = ^ n 


b{i) 


E 


1 


/c=i + + 

The asymptotic fitness distribution is 

k-l 


C 


Pk = 


n 


i + l 


A* + (1 + (3){k + 1 ) A* + (1 + (3){i + 1) 


= C 


, r(A: + l) 


(t + 3+5^) (l + «‘ 


t-t+ral 1 


(l + /9)‘ 




A mechanism where the vertices die at a rate proportional to their htness 
hence destroys the power law distribution obtained without vertex deaths 
regardless of how small the proportionality constant is and instead produces 
a (super)exponentially decaying htness distribution. 

We remark that b{i) = (z + 1)'^ and d{i) = I3{i + 1)'^ gives rise to the same 
exponential decay for any 7 > 0. Indeed, 


_ \ _TT_ + 1)'^ _ 

A* + (1 + (3){k + 1)7 A* + (1 + (3){i + 1)7 

^ e“ l°g(i+(i+/3Ki+i)'/) 

(1 + /3)" 
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where the sum in the exponent grows like k} for 7 < 1 and is convergent 
for 7 > 1 . 


Example 3. Next take b{i) = i + 1 and d{i) = + 1)^ for /3 < 1 and 

rj € ( 0 , 1 ). Also in this case p{X) < 00 for all A and, by comparison with the 
previous example, (3 < 1 guarantees that the process is supercritical. For the 
fitness distribution we get that 



Vertex deaths determined by a sublinear function of the fitness hence also 
destroy the power law distribution obtained without vertex death, bringing it 
down to a stretched exponential distribution. In contrast to the case when 
the vertices die at random however, the tail behavior of the distribution is 
changed by the vertex deaths. 

Again we remark that b{i) = {i + 1)'^ and d{i) = f3{i + 1)*^ for arbitrary 
7 > r/ > 0 with 7 —r/ € ( 0 , 1 ) gives rise to the same type of htness distribution: 
a similar derivation as above implies that logp^ ~ 

Example 4. Finally, consider a situation with b{i) = a and d{i) = (f + 1)~^, 
that is, new vertices attach randomly to existing vertices and vertex deaths 
occur inversely proportional to the htness. Analogously to ([9]) we get that 



and using (l 8 |) it follows that 



Hence the process is supercritical if and only if a > 1/2. The Laplace trans¬ 
form is 



where again F is the hyper geometric function, and the equation p{\*) = 1 is 
easily solved numerically for a given value of a. As for the htness distribution. 
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we have 


Pk = 


k-l 

:rn 


a 


A* + a + ^ A* + a + ^ 


a 


fc+i 


r(A; + 2) 


a\\* + a) r(fc + 2 + ^) 


a 


A* + a 


k ^*+“, 


that is, the fitness distribution decays (super)exponentially at rate a/(A* +a). 
The fact that vertices that are not successful in acquiring edges dies at a higher 
rate is hence not sufficient in itself to explain the power law distributions for 
the fitnesses. 


5 The in-degrees 


As mentioned in Section 1.1, the fitnesses (that is, the accumulated in-degrees) 
are not seldom the primary object of interest when studying a network incor¬ 
porating vertex death. However, there are also situations when the momenta- 
neous in-degrees, henceforth referred to as just the in-degrees, are the relevant 
quantities. This is presumably the case for instance when modeling the World 
Wide Web. In this section we demonstrate how a formula for the asymptotic 
in-degree distribution follows from the same method as for the fitnesses and 
present simulation results that show that the in-degree distribution behaves 
in a similar way as the fitness distribution. 

Consider a vertex population that evolves in the same way as described 
above, with the difference that the birth and death rates are based on the 
in-degrees rather than the accumulated in-degrees. More specifically, a vertex 
that is alive at time t and has i alive children at time t gives birth to a new 
vertex at rate b{i) and dies at rate d{i). Let ^ be a random point process 
on [0, oo) with points (representing the births) generated according to these 
rules. We use the same notation as before equipped with a wave-hat for 
quantitites related to omitting the original hat for the Laplace transform of 
the reproduction function. The Laplace transform is given by 


P(A) = '^b{k) 

k=0 



k n L > t] dt, 


( 10 ) 


and A* is defined by p{X*) = 1. Furthermore, if b{i) and d{i) satisfies the 
analogues of (C.l) and (C.2), it follows from Theorem 12.11 that 



t^OO 


/o°° 6*^*^ (ii't) = kn L > dt 

/o“P(i>t)dt 


( 11 ) 


13 










that is, the fraction of vertices with in-degree k converges almost surely and 
the limit is given by the right-hand of m- Unfortunately it does not seem 
possible to derive a more explicit expression for the limit in this case. The 
event that a vertex of age t has k alive children occurs if the vertex has given 
birth to n children for some n > k and n — A: of these children have died before 
the vertex has reached the age t. To find an expression for the probability of 
this event that leads to a tractable formula for the integral in (jlOh and in the 
numerator of (fTT]l does not seem possible. 

To get a grip of the asymptotic in-degree distribution, we turn to simu¬ 
lations. These indicate that the qualitative behavior of the in-degrees is the 
same as for the fitnesses. Figure 1 shows simulated in-degrees and fitnesses 
plotted on a log-log scale for b{i) = i -\-l and d{i) = 0.5. In Example 1 it was 
shown that the asymptotic fitnesses in this case follow a power-law distribu¬ 
tion with exponent 3, and Figure 1 reveals that the in-degrees seem to follow 
the same type of distribution, although the exponent appears to be slightly 
larger. Figure 2 shows simulated in-degrees and htnesses for h{i) = i -\-l and 
d{i) = 0.5(i -I- 1), with a log-scale on the y-axis. The asymptotic fitness dis¬ 
tribution was shown in Example 2 to have an exponentially decaying tail in 
this case, and the same seems to be true for the in-degrees. Finally, Figure 3 
shows simulated in-degrees and fitnesses on a log-log scale for h(i) = i-\-l and 
d{i) = 0.5(i -|- The plot reveals a slightly curved behavior for both the 

in-degrees and the htnesses which is signihcative for a stretched exponential 
distribution and in agreement with the analytical result for the htness dis¬ 
tribution. Note that, in all three cases, the tail of the in-degree distribution 
seems to decay faster than the tail of the htness distribution. All simulations 
have been run until the number of alive vertices equals 100,000. 

6 Further work 

There are a number of questions about the studied model that deserve further 
investigation. What about the global properties of the obtained structure for 
instance? If the ghost links are included one gets a random tree, analyzed in 
Rudas and Toth m, and excluding the ghost links gives rise to a forest. Will 
there be an inhnite component in this forest? How does the answer depend on 
b{i) and d{i)7 Another issue to look into is the time evolution of the process. 
For models without vertex death it is known that the maximal degree grows 
polynomially of the same order as the first vertex; see Mori m- Is this the 
case also in the current setting for the maximal htness? Recently Dehrich 
and Morters m have shown for a model with attachment function k"^ and 
no vertex death that there is a phase transition at 7 = 1 / 2 : if 7 >l /2 there 
is a persistent hub that will be of maximal degree at all but at most hnitely 
many times, while, if 7 < 1/2, there is no persistent hub. Is there a similar 
phenomenon in a model with vertex death? 
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Figure 1: Simulation with b{i) = z + 1 and d{i) = 0.5. The figure shows the tail 
probabilities of the fitness distribution (stars) and the in-degree distribution 
(circles). 
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Figure 2: Simulation with b{i) = i + 1 and d{i) = 0.5(z + 1). The figure 
shows the tail probabilities of the htness distribution (stars) and the in-degree 
distribution (circles). 
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Figure 3: Simulation with b{i) = i + 1 and d{i) = 0.5(z + 1). The figure 
shows the tail probabilities of the htness distribution (stars) and the in-degree 
distribution (circles). 
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